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CHAOTIC DISTRIBUTIONS FOR RELATIVISTIC PARTICLES 


D. MUSTAFA AND B. WENNBERG 


Abstract. We study a modified Kac model where the classical kinetic 
energy is replaced by an arbitrary energy function cp(v), celR. The aim 
of this paper is to show that the uniform density with respect to the micro- 
canonical measure is Ce~ Zo ^h)- c haotic, C,zq £ R+. The kinetic energy for 
relativistic particles is a special case. A generalization to the case v 6 R rf 
which involves conservation momentum is also formally discussed. 


1. Introduction 

In 1956, Mark Kac published a paper 0, in which he answered some 
fundamental questions concerning the derivation of the spatially homogen¬ 
eous Boltzmann equation. Kac considered a stochastic model consisting of 
N identical particles, and obtained an equation like the spatially homoge¬ 
neous Boltzmann equation as a mean-field limit when the number of par¬ 
ticles tends to infinity. The key ingredient in his paper was the notion of 
chaos, which goes back to Boltzmann's stosszahlansatz. Loosely speaking, 
chaos is related to asymptotic independence. We give a precise mathemat¬ 
ical definition later. He showed that if the probability distribution of the 
initial state of the particles is chaotic then this property is propagated in 
time, i.e., the probability distribution at time t > 0 is also chaotic. This is 
referred to as propagation of chaos. 

Kac's model is a jump process on S N “ 1 = {(zq,..., zqy) £ 1R N | v 2 + • • • + 
v 2 n = N}. This corresponds to an N-particle system where a particle is rep¬ 
resented by its velocity v G 1R, and the kinetic energy of a particle is given 
by v 2 . The total energy is conserved in this process, which is ergodic and 
hence has a unique invariant measure. In the original model, the station¬ 
ary distribution is the uniform measure on S N ~ 1 . It is easy to prove that 
this uniform measure is (2rc) ^ 2 e v //2 -chaotic according to the definition 
given below, and Kac also provides a method of constructing other families 
of chaotic measures on S v 1 which may serve as suitable initial distribu¬ 
tions for the process. 

In this paper we consider a generalization of Kac's model, where the en¬ 
ergy of a particle is given by a large class of functions (p(v), which includes 
(p(v) = v 2 . The motivation for studying this problem comes from the theory 
of special relativity, where v represents the momentum of a particle rather 
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than its velocity, and where the kinetic energy is given by 

(p(v) = \J\ + v 1 — 1. 

The kinetic theory for relativistic particles is far less developed than for 
classical particles, but there is a growing interest in this kind of problems 
(see below for references), and this paper is a contribution to this work. 

The first result of this paper is a proof that the uniform distribution on the 
manifold represented by the total kinetic energy of the particles (depends 
on (p ) with respect to the microcanonical measure is chaotic by using the 
approach of Kac. We also consider a jump process on this manifold, which 
is a direct generalization of Kac's model, and prove that propagation of 
chaos holds for this process. We also consider a generalization to a model 
where the particles have momenta » £ E‘l 

One of the most influential equations in the mathematical kinetic theory 
of gases is the Boltzmann equation, which describes the time evolution of 
the density of a single particle in a gas consisting of a large number of par¬ 
ticles. 

The Boltzmann equation is 

Vx/(x,u,f) = Q(f,f)(x,v,t), v £ R 3 ,f > 0, 

\/M)=/o, 

where in the classical case v = v, and in the relativistic case v = v/ ^Jl + \v\ 2 . 
The collision operator on the right hand side is a quadratic operator acting 
only in the velocity variables. In the classical case it is given by 

Q(f,f) = J r 3 J s2 -/(>,f)/(T*,f)) B(v - v*,cr)dv*dcr. 

The function B is called the collision kernel and is derived from the physics 
of the model. The corresponding collision term is considerably more com¬ 
plicated in the relativistic case, which is thoroughly investigated in lfl5Tl . for 
example. In both cases, (v, u*) and (v 1 , v[) represent the velocities of a pair 
of particles before and after an elastic collision, and satisfy the conservation 
of momentum, 

v + u* = v’ + z/*, (1.2) 

and energy 

<p(\v\) + <p(\v*\) = <p(\v'\) + <p(\v'* |). (1.3) 

The derivation of the Boltzmann equation from a classical (Newtonian) 
many particle system is still an important research topic. The classical result 
in this direction is by Lanford ITOl . however, this result is valid only over a 
fraction of the mean free time between collisions. Up to this date the result 
by Lanford has only been improved upon only in details, see e.g., 1121 and 
the recent book 0. For the study of the relativistic Boltzmann equation, see 
e.g., El. 
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To obtain an equation like (11.11) as a mean-field limit, Kac considered the 
master equation, which describes the evolution of the probability density on 
state space. Assuming that velocities of the particles is distributed accord¬ 
ing to a probability density Fq initially, its time evolution is given by the 
following linear master equation 

f . .,VN,t) = K.F N (vi,. . . ,1’N,t), ^ ^ 

\f n (vi,...,v n ,0) = F 0 (v lr ...,v N ). 

where K. is the generator of a Markov jump process, 

fcF N (V\,.. .,v N ) = N[Q - I]F n (v i,.. .,v N ), 

with I being the identity operator and 

2 f n d6 

QF N {v v n ) = 1)^.1 _ n F N{vi, ... rV'iiO),... ,v'j{0),... ,v N ) — . 

Hence /C can be written as a sum of generators acting on only two variables, 
Vj and Vj, say, which corresponds to binary collisions in a real gas. The pair 
of velocities v\[0),v'-(0) is the outcome of the velocities v u Vj undergoing a 
collision, and in the classical case, as proposed by Kac they are given by 

v’^0) = Vi cos 9 + Vj sin 6, v '(0) = — u,sin0 + vjcosO. (1.5) 

The energy of a pair of particles is always conserved in a collision, i.e., 

v'(e) 2 + v'(6) 2 = vj + vj. 

In a realistic model the momentum should also be conserved, but since the 
particles have one-dimensional velocities, imposing conservation of mo¬ 
mentum would lead to either the particles keeping the velocities during 
a collision or exchanging velocities. Momentum conservation is therefore 
sacrificed in this case. 

The fact that particles are assumed to be indistinguishable corresponds 
to the initial distribution Fo being symmetric with respect to permutations, 
and this property is preserved by 1C, so that F ( \; is symmetric for all t > 0. 

In order to obtain a mean-field limit of the master equation (11.41) as N —> 
oo, Kac introduced the notion of chaos (in 0, he referred to it as the Boltz¬ 
mann property). 

Definition 1.1. Let / be probability density on K. For each N £ N, let 
F/v be a probability density on IR :V with respect to a measure m IN} . The 
sequence {Fn}n eN of probability densities on 1R ;V is said to be /—chaotic if 
the following two conditions are satisfied: 

(1) Each Fn is a symmetric function ofvi,..., Vm- 

(2) For each fixed k £ N, the k- th marginal f^(v \,..., ly) of Fn con¬ 
verges to nil f{ v i)> as N —t oo, where f(v) = limN-^oo/ 1 N (f). The 
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convergence is to be taken in the weak sense, that is, if q>{v \,..., v k ) 
is a bounded continuous function of k variables, G IR, then 


lim / cp(v 1 ,...,v k )F N (v 1 ...,v N )dm {N) 

N-foo JrN 


= / 

J id 


■,v k ) n/( Vj)dv i... 

z=l 


( 1 . 6 ) 


More generally, this can be formulated for probability densities on prod¬ 
ucts E n where E may be an arbitrary Polish space. The setting of this paper 
is also slightly different from the Definition 1.1 because the dynamics take 
place on a submanifold of R N . 

A chaotic family of probability densities on S N_1 (VN) is the following: 


Example 1.2. It is a well known fact that the surface area of the sphere 
S N " 1 (V'I) in R N is given by 


|s n - 1 (Ve)| = 


2ttteV 


Let 


Fn(v = 


zN-2 


E — v\ 


dv i. 


|SN-1(^N)| 

be the symmetric uniform density on S N “ 1 (v / N) with respect to surface 
measure a i N K Let cp be continuous function on R. It is easy to see that 

yjE 

( cp(v 1 )dcr {E ') = [ (p (zq) 

J SN-i(v'E) J-Ve 

Replacing E with N we have 

/sn-i(Vn) <p( y i )dr {N) 

lim--- - -—- 

n -»» |s n ~ 1 (Vn)\ 

L(f) 


= lim 


r-v/N 


N—3 


<p(Ul) 1 - 


N 


dm 


1 F 00 

Taking <p to be a function on R /c , the same calculations show that, the family 
of uniform densities on S :V 1 ( \f~N ), with respect to the surface measure, is 

1 _ 1 
-4=e 2 -chaotic. 

v27T 

Kac showed by using a combinatorial argument that the master equation 
(11.41) propagates chaos, that is, if {F^(v i, ... ,vn, 0)}nsn is/ q— chaotic, then 
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the solutions to (11.41) . {F N (zq, ..., Uy, t) } N€]N is f(v, t)— chaotic. The density 
f(v,t ) satisfies the Boltzmann-Kac equation 


i) r r n dO 

gj/M) = 2 j ]R J_ 7 i \f(v(e),t)f(w( 0 ),t) - f(v,t)f{w,t)\ —dw, 

f{v, o) =fo, 


with v(6), w(0) given by (11.51) . 

Kac described also a method to construct other chaotic probability den¬ 
sities on S N - X Let 


Fn{v i,.. .,Vn) = 


nfuM) 

z(Vn) ' 


where 


r N 

Z(VE)= n f(?i)dv {E \ 

Js n -'(Ve) i = i 


( 1 . 8 ) 


(1.9) 


with a 11 ; being the surface measure on S n ~ 1 (VE). Kac showed that, un¬ 
der some conditions on the function /, the family of probability densities 
{FnInen is h(v) —chaotic, where 


h(v) = 


-z 0 v 


f(v) 


/r e ~ z ° 1 ’ 2 f (v)dv' 


( 1 . 10 ) 


and Zo is a positive constant. Note that the case f(v) = 1 corresponds to the 
previous example. 

The purpose of the present article is to study chaotic probability densities 
on other surfaces than S ;V 1 (VN). Define 


cp(vi) = Ej, (1.11) 

where (p(v ) represents the energy of a particle with velocity v, and 

cp G C 1 (R,R + ) / (p is even and convex, 

<p( 0) = 0. (1.12) 

We show that the uniform density on CW 1 ( \f~N) with respect to the mi- 
crocanonical measure is Ce~ z ° < f > '- >— chaotic, where C is a normalisation con¬ 
stant and Zo > 0 is the unique solution to a specific equation. 

Chaotic probability densities on spheres have been considered by many 
authors. In a paper by Carlen, Carvalho, Le Roux, Loss and Villani, |l|, the 
authors show the following theorem, which differs from Kac result both in 
the exact definition of chaos (the notion of entropic chaos is introduced), 
and in the method of proof: 


N 


n N 1 (Ve) = < | 


i =1 
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Theorem 1.3. Let f be a probability density on R such that f G L? (R) for 
some p > 1, f ]R v 2 f(v)dv = 1 arid f R v 4 f(v)dv < oo. Let a family of densi¬ 
ties {FjvjiveiN/ be defined with F ;v as in il.8t with now being the normalized 
surface measure on S n ~ 1 (i/N). Then {F]v}n6N is f-chaotic. 

Sznitman lfl3ll constructed chaotic families of measures using a some¬ 
what different approach: Let hi ,..., h\< be i.i.d random variables with law 
y(dh) = f(h)dh, where h G R d . Assume that / is differentiable and satisfies 
the following condition 


f if(h) + \Vf(h)\)e m dh 
J R rf 


< 00 , 


(1.13) 


where A G R. Then, the conditional distribution u N of (hi,..., /tv) subject 
to the constraint /;i+ ’ w +/ ' N = a is Y—chaotic, where 



K-h 


p(dh), 


(1.14) 


with A„ determined by the equation J hdY(h) = a and is a normalisation 
constant. Note that, the choice of A for which condition (11.131) is satisfied 
depends on a. Within the framework of this paper we think of the random 
variable It G R as the energy of a particle, in this case h = <p(v) >0. 

The organization of this paper is as follows. In section 2 we define the 
microcanonical measure on Q N 1 ( \f~N) and show that the uniform density 
on _ 

n N -\VN) with respect to this measure is chaotic. In Section 3, we intro¬ 
duce a modified Kac model with particles having the energy <p(v) and the 
corresponding master equation. We also show that this master equation 
propagates chaos and obtain the limiting equation. In Section 4, we dis¬ 
cuss how to generalize the results in Section 2 to the case v G R d , d > 1, in 
which case the momentum is also conserved. An Appendix is also included 
to introduce the methods that are used in the paper. 


2. Chaotic measures 

This section is devoted to showing that the family of uniform probabil¬ 
ity densities on f) v 1 ( x^N ) with respect to the microcanonical measure is 
chaotic. By Example ll.21 another approach would be show that the uniform 
density on Q N_1 (\/N) with respect to surface measure is chaotic. However, 
since we consider a particle system where the total energy is conserved it 
is natural to use the microcanonical measure on Q v l (\^N). We show that 
the family of uniform probability densities on C1 n ^ 1 (i/N) with respect to 
the microcanonical measure is Ce~ z ° ( l ’^-chaotic where Zq > 0 is the solution 
to a specific equation and C is a normalisation constant. 





CHAOTIC DISTRIBUTIONS FOR RELATIVISTIC PARTICLES 


7 


The microcanonical measure on Q ;V 1 (VN) is defined as follows: Let 

N 

H(vx,...,v n ) = Y^(p{vi). (2.1) 

i=i 


Definition 2.1. Provided that | VH| ^ 0, the microcanonical measure, i] iE) , 
on Q N “ 1 (v / E), is defined by 


V 


(E) = 


\VH\' 


( 2 . 2 ) 


where <tq is the surface measure induced by the Euclidean measure in IR N 
on Q n ~ 1 (x/~E). 


The microcanonical measure arises naturally in physics from the assump¬ 
tion of equal probability, meaning essentially that all microscopic particle 
configurations corresponding to the same energy are equally probable, see 

e.g. ini or mg. 

A more geometrical approach is given by the coarea formula, as explained 
e.g. in |6): If <f> : IR" —> R m is Lipschitz continuous and n > rn, then for each 
measurable set A C IR", 

[ J<P{x)dx= [ U n ~ m {A^^ l (y))dy, (2.3) 

Ja Jw n 


where JO is the Jacobian of O, as defined in J5J, and J{ n ~ m denotes the 
n — m-dimensional Hausdorff measure. In our setting this can be rephrased 
as (see Theorem 3.13 in f6|): if H : R" —» R is Lipschitz continuous and 
ess inf| VH| > 0, and g : R n — > R is integrable, then for almost all E, 

Je L £MVhm>e) * = f [H=E} ■ (Z4) 

where x is the indicator function of a set. Because H(v\, ...,Vn) = cj)(v i) + 
■ ■ ■ <P(vn) with ([> convex and diffentiable, the set {H = E} is also convex 
and regular. Then 


/ g(v i,..., v n )S(H(v i, ...,v N )~ E)dv i... dv N 
J]r n 


R2nw-i(yy) 




ln*-\VE) 8 ]VH\- 




This shows that 

5(H(v lt ... / v N )-E) = ^ f 

Another concept that is relevant in this context is that of disintegration 
of a measure, which is a measure theoretic approach to almost the same 
problem, and a means of approaching conditional probabilities, see e.g. |5). 
In our setting, the disintegration theorem states that given a measure dy on 
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R n there are a measure v on R + and a family of measures ^ on the level 
surfaces {(vi,...,Vn) \ T,(p(vj) = E}, such that 

[ N f(y)M u )=[ f JWmWviE ), ( 2 . 5 ) 

J r n Jo J{z<p{vj)=E} 

Remark 1. On S n- 1 (\/E), we have that |VH| = 2\jE. This implies that 
the microcanonical measure is up to a constant factor equal to the surface 
measure on S n_1 (\/E)- 

Using the equality H(v . .., V\/) = E, we can express, at least locally, the 
variable Vn as function of V\, ... , v^-i- 

v N = U(v i,...,y N _i). 

By this parametrization, the surface O v 1 (-\/E) can be represented as the 
graph of U : R N ~ X —> R. The surface measure an on Q n_ 1 (\/E) is now 
given by 

dan = + | Vlt| 2 dv\... dv^~ i- 

By the implicit function theorem it follows that 

dH 


= k = -l . N-l. 


Thus 


dan 

Wm 


dH 


dv N 


dv\... dv^-i. 


To carry out integration on Cl N 1 (\/E) with respect to the microcanonical 
measure ;/ (£) we use the last equality: 


/nw-pv/E)' 


= E 




g(vev N ) 


dH 


dv N 

€ 


dv i... dvN-i, 


where 


vn = (P 1 ( E ~ E / 


and (p 1 (v) is the of cp(v ), v > 0. Moreover, 


dH 


dv N 


dH 


N-l 


dv N 


v 1 ,...,v N - 1 ,e<p ME- E <P( V i) 
v V i= 1 


The uniform density F(v i,..., z?n) on Q N 1 (VN) with respect to the mi¬ 
crocanonical measure if Nj is given by 

Fn(v ' . Vn) = (Z6) 
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where 


^(s/E) = [ dt,(E). 

H J nN-i(VE) 


To show that the uniform density on C1 N with respect to the micro- 

canonical measure i] lN} is Ce~ z °^^ v )-c haotic we follow Kac (9], and start by 
determining the asymptotic behaviour of 

-dv i... dujsi -1 


= £ 

'E : 1 E 


dH 


dv N 

e 


with E = N for large N. Since 
dH 


dv N 


= \(p'{v N )\ and = ±</> ME- <K^) , 


N—1 

E 

i=l 


we have 
ZA\ 


[ 

l 

'eMMM)<e 



... dvj\i— 


To write Z^(v / £) as an integral over the sphere S N 1 (-\/E), we make the 
change of variables y\ = <p(vi ) with respect to sign of v u i = 1,..., N — 1. 
This leads to 


z^Ve) = 




JrFn 1 


N—l 


\yi 


yf<E 




SI W(.<P~ l {y}))\ 


Let 


/(y) := u 77 


|y| 


dyi... di/N-i- 

(2.7) 

( 2 . 8 ) 


l^'(r 1 (y 2 ))r 

The integrand in (12.711 is almost a product of N copies of /(y). Multiply 


and divide the integrand by |yM = \J E — Recall the following 

formula for integration over a sphere 


/ g{yi'■ ■ ■ fVN)dcr (E ^ 

Js n -He) 


> f / \ Edy\... dy^—i 

-SM**** . 


We now get 


9N-I N 


r(E) 


(2.9) 


Having Z<p given by (12.911 is convenient in sense that, in 13, Kac deter¬ 
mined the asymptotic behaviour of Z^(v / N) for large N by using the saddle 
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point method (see e.g. 03). For completeness, we present each step of the re¬ 
sult with rigorous justification with f{y) given by (12.81) . A short description 
of the saddle point method is given in the Appendix. 

We start by computing the Laplace transform of E H* Z ( p(y/E). The 
Laplace transform of Z^{\[E') is defined provided that Z,p (V£) grows at 
most exponentially. Since the behaviour of Z l p(y/E) depends on the func¬ 
tion f(y) defined by (12.81) we assume that cp(y) is such that 

f(y) < Ke by2 , ( 2 - 10 ) 

for some K > 0 and b > 0. This condition ensures that Z f p(\/£) grows at 
most exponentially. 

Taking the Laplace transform of Z^^y/E), making the change of variable 
r = y/E, we have, for iv 6 C where 9t(w) > b 

nOO _ nOO 

e~ zvE Z (p {VE)dE = 2 j e~ wr rZ^{r)dr. 

Using (12.91) . the last equality equals 

2 N y\--rf(y)d y y . 

From condition (12.1011 and since 3J(zy) > b it follows that 

O(za) := r e-^ 2 f(y)dy (2.11) 

J —CO 


is an analytic function of w for Sft(ic) > b. By applying the inverse of the 
Laplace transform, we get 

2 N rj+ioo / poo 2 \ N 

Z${Ve) = — 1^ e z yj J zy f(y)dyj dz, 

where 7 > b is such that the line 7 = jft(z) lies in the half-plane where 0 ( 2 ) 
is analytic. Replacing E with N, we get 


z^Vn) 



N 

dz. 


By comparing with the saddle point integral d5.il) in the Appendix, we set 


q(z) = l and S(z) = z + log <3>(z), (2.12) 

which is well defined because <t>(z) is nonzero on the line 7 = R(z). We 
now have 

oN-1 /* / y+ioo 

Z<p(Vn) = -- / e NS ^dz. (2.13) 

7X1 Jy—ioo 

The asymptotic behaviour of Z l p(y/N) for large N is determined by the sad¬ 
dle points of S(z). The next lemma concerns the saddle points of S(z). 
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Lemma 2.2. Assume that there exists a 7 > 0 such that 

[ e~ iy2 fWy < 00 . 

J —00 


Moreover, assume that 


and 


/ OO 

(i - y 2 )f(y) d v < o, 

-00 

l , (i - y 2 )f(y) d y > o. 

J\v\< 1 


(2.14) 

(2.15) 

(2.16) 


Lef S(z) be given by (I2.12D with <f>(z) fey (12. IIP . For 9ft(z) > 7 , the function S(z) 
is analytic and there exists a unique saddle point zq to S(z) such that zo is reaZ, 
z 0 > 7 and S"(zq) > 0. Moreover 


z#{Vn) 


2 N-l e Nz 0 


(/» e zoy2 f(v) d yj ■ 


VNS"(z 0 ) 

Proof. Note that, for z = £ + zy, for all £ 

arg max |e?+"/<t>(£ + iy)| = { 0 }. 
v 

Hence, we only need to find saddle points on the real line. 

Claim V. The function S(z) has a unique saddle point Zq where Zq > 7 . 
Proof. For 'ft(z) > 7 , we have 

S-ooV 1 ^ 2 f(y) d y 


(2.17) 


S'(z) = 1 - 


(2.18) 


J-^Mdy • 

For z = f + ZO where £ > 7 , the equation S'(£) = 0 is equivalent to 

/ °° 9 

e -c*/ (1 _ y 2 )f( y )dy = 0 . 

-OO 

Multiplying the last equality by ef we see that, S'(f ) = 0 is equivalent to 

/ CO 

-OO 


Let 


N(y 2 -i)n _ iy 2 


(i - y 2 )f(y) d y- 


Since 


r 

= J 

/ CO 

e -f(y -D( 1 _ y 2 ) 2 /(j/)rfy>0/ 

-OO 

it follows that A(£) is an increasing function of £. Moreover, by (12.151) we 
have 


Note that 


lim [ 

£->°° Jhi 


21 ( 0 ) < 0 . 


N(y 2 -i)n_i/2 


>1 


(1 - y 2 )f(y) d y = o. 
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Hence 



^-Dfl-/)/(y)dy 



e -«J J -l) ( l_ y 2 )/(y) rf ¥ 


= f <'■" 1 - ihflytdy. 


The last integral goes to infinity by (12.161) as £ —> oo. Hence, there exists 
unique Zp > 7 such that S'(zp) = 0. 


Claim 2: The second derivative of S(z) at Zq is positive. 


Proof. We have 

_ *"(*) <^(z) 2 

1 j O(z) <3>(z) 2 ' 

For z = Zo, using the Jensen inequality, we get 

<p"(zo) = /-«, y 4e ~ Zolj2 f(y) d y ( IZ>y 2e ~ Zoy2 f(y) d v \ 2 = ^(zp ) 2 
°( z o) J-co e ~ zoy2 f(y) d y v /“oo e~ Zoy2 f(y) d y ) °( z o) 2 ' 

This proves the claim. □ 


We now turn to the proof of (12.17b . We can write 

O(z) = / e~ zy2 f(y)dy = f e~ zy2 (/(y) + f(-y))dy. 

J — oo JO 

By a change of variables, the last integral equals 

- 2 V /(>/y) + /(- v^). 


/ 

do 




-dy. 


Let z = £ + ///, with £ > 7 . The last integral can be written as 


/ 

./ —I 


2v/y 


■l y >o e " /y dy. 


The last integral is the Fourier transform of the function at the point //, 
where 


&e(y) 


£ J(Vy) + f(-Vy) 

2 


ly>0- 


For £ = Zo, it follows that Og E L 1 (R), and | ( 77 ) | < More¬ 

over, By the Riemann Lebesgue lemma, it follows that 


|J 7 (Og)(y)| —> 0 when \rj\ —» 00 . 

Let C + t and C 7 be the curves in the complex plane given [7 + zT, zo + zT] 
and [7 — zT,zo — zT], respectively. 

We have 



= [~° (T) N d£ 

Jy 

rZ 0 

< / e N ^(^)(T)| N ^ 

■n 


□ CD 
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< \zo — j\ max e N ^\ J r (0^)(T)| N —» 0, as|T| —> oo. 

fe[7,zo] 


By Cauchy's theorem, we can deform the contour in (12.1311 to a contour 
passing through the saddle point zq. Hence 

oN-l rZofyioo 

ZJVn) = / e NS<yZ ^dz. 

TCI Jz 0 -ioo 

By the saddle point method, we obtain 


z^Vn) = 




2^-1 
n i 



1 + 0 



e NS(z 0 ) 


2 n £ Nz 0 


(sy*fw* 


(2.19) 


This finishes the proof of the lemma. 


□ 


The main theorem of this section is: 


Theorem 2.3. Let f(y) defined by (12.8I > satisfy the conditions (12.101 ). (12.151 ) and 
(I2.16D . Then, the family of uniform densities on Q n ~ 1 (-\/N) with respect to the 
microcanonical measure is Ce -20 ^ 1 ^ — chaotic. The positive constant z$ is the 
unique real solution to the following equation 


/ OO 

(1 — cp(v))e~ Zo< l ,l ' v ' > dv = 0 , 

-oo 


and C is a normalisation constant given by 


f°° e~ z o<P( v )dv' 

J —OO 

Proof. Let cp be a bounded continuous function on IR^. For each fixed k G IN, 
a modification of Lemma fZ2l shows that, for large N, 


V / 


cp(vV k )- 


dH 

dv N 


-dv\ ... 


m k 


cp((p l (y\),...,(p 1 {y 1 k ))Y[e z °y 1 if{y i )dy 1 ...dy k 


i=i 


2 N e (N-k)z 0 

VNS"(z 0 ) 
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Using Lemma 12.21 again and making the change of variable y 2 = <p(vi), 
where dvi = f(yi)dyi, i = \,... ,k leads to 


lim 

N —>oo 


InN-i {V N ) (p^ir---,Vk)dii iN) 

Z<p(\/N) 

f. ?(0 _1 (y 1 )/ • • ■ '<P~ l {yl)) Y[e~ Zoy ‘f(,yi)dyi ■ ■ ■ dy k 

jR k 


/Ri i =1 


Yl e Zoy ‘f(yi)dy-\ ■■■dy k 

i 

k 

1 

i =1 


[ (p(v-i,...,v k )Y[e dv 1 ... dv k 
J R* 7-T 


, k 

/ • • • dv k 

J ]R* + = Y 


This is what we wanted to prove. 


□ 


We now consider two examples where Theorem l2. 31 applies : 


Example 2.4. In the classical Kac model (p(v) = v 2 . We get that f(y) = ^ 
and thus, the conditions (12.141) . (12.151) and (12.161) are fulfilled. To determine 
Z() we need solve the equation 



z ° v2 dv = 0 . 


Direct calculations show that Zq = \ is the unique real solution. Therefore 
the uniform density on (Q :V = S v 1 (v^/V J) with respect to the mi- 

_l v 2 

crocanonical measure is e , — chaotic. This has already been discussed in 


Example 11.21 Recall that, on S N 1 (), the microcanonical measure up to 
a constant is equal to the surface measure 


Example 2.5. For a relativistic energy function (p(v) = \Jv 2 + 1 — 1, it fol¬ 
lows that 

(y 2 + 

vf+FT 

It is easy to check that the conditions (12.141) . (12.151) and (12.161) are satisfied. 
To find Zq we solve the equation 




/ OO _ 

(1 - (vV + 1 

-OO 


1 )) e -2 0 (x^ur-i)^ = a 


By using numerical integration, we get that Zo ~ 0.734641. Hence, the fam¬ 
ily of uniform densities on Q n_1 (-\/N) with (p(v) = Vv 2 + 1 — 1 in dl.lll) is 
Ce -zo(V?+T-i)_ chaotic, where C ~ 4.082. 
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We end this section by a comparison between our result and the result 
of Sznitman.The setting there is slightly different, and as commented by 
Sznitman it is not directly applicable for the case h = v 2 or the more gen¬ 
eral situation here. To see the relation between two results, assume that an 
individual velocity is distributed according the Gibbs measure /Z, 

where Z is normalization constant and /3 > 0. Then the random variable 
h = (p(v) that is the energy of a particle, has a distribution y(dh) = f(h)dh, 
where 

e ~ph 

= WlrHM' (2,20) 

To apply Sznitman's result, we need to show that / satisfies condition 
HI.1311 . This is a strong integrability condition on / and depends on the 
choice of cp . Unfortunately, even for the most classical case, the one studied 
by Kac, where h = v 2 , condition (11.131) is not fulfilled since f'(h)e~P h is 
not integrable at h = 0. It is a technical condition needed to control the 
Fourier transform of /, and it could probably be relaxed, but it means that 
Sznitman's result cannot be directly applied to our case. 

3. Particle dynamics and master equation for general energy 

FUNCTIONS 

In this section we follow |2] and f9] to introduce dynamics between par¬ 
ticles having the energy given by the function (p and obtain the correspond¬ 
ing master equation. By similar arguments as in 03, we will see that the 
master equation propagates chaos. 

To introduce dynamics between the particles, let the master vector 
V = (vi,... ,Vn) G O n_ 1 (v / N). The master vector V makes a jump on 
O N 1 ( \/~N) according the following steps: 

(1) Pick a pair ( i,j ), i < j according to the uniform distribution 

p.. = _ - _ 

11 N(N - 1) ‘ 

(2) The pair of velocities ( Vi,Vj ) satisfies 

cp(vi) + cp(vj ) = h, h > 0 . 

Let 

Vi = sign(u) \!<p{vi) and \j } = sign (v,) ^<p{vj). 

Then (y,-, ijj) is a point on the circle, and as in the original Kac model, 
the collision may be performed there. 

Pick an angle 6 uniformly on (0,2zr] and let 

y■ ( 0 ) =\ji cos 6 + ijj sin 6 and y' ( 0 ) = —y, sin 9 + y l cos 6. 

The pair ( Vi,Vj ) is transformed on Q 1 (v / /i) to (v' i (9),v'j(6)) according 
to 



v'i(0) =sign(y'(0))<£ 1 (y'(0) 2 ). 
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where 


Note that 


^■( 0 ) = si § n (y/( 0 ))^ 1 ' 


sign (v) 


1, if v > 0 

— 1, if v < 0. 


(3.1) 


cp(v'i(d)) + <p(v'j(e)) = <p(Vi) + <p(Vj). 

(3) Update the master vector V and denote the new master vector by 
Tj r j(9)V. By step 2, it follows that T,y (0)V 6 Q ;V1 (\/7s7). Repeat 
step 1, 2 and 3. 

This is only a generalization of the dynamics in (9) where <p(v) = v 2 . 

The steps above describe a random walk on Q n_1 (\/N). As in 0, we 
define its Markov transition operator Q,p. If V/ c is the state of the particles 
after the k-th step of the walk, for a continuous function cp on Q ;V 
the operator Q ( p is defined by 

QM y) = e W(v k+ i)|v3k = y]. 

Writing out the expectation above, we get 

= JW^Yi^r <3 ' 2) 


If Fj, is probability density of V/- with respect to the microcanonical measure 
yt'U on n N ~\VN ), we have 

/ nN _ )^ + !^ (N) = E[^(V t+1 )] = E [E[(p(y fc+ 1 )|Vfc]] 

= /' ^ QftpFkdrjW. 

Jci n - 1 (Vn) 

By definition, the microcanonical measure is invariant under the transfor¬ 
mation V —> T,y (0)V. It follows that Q^, is self adjoint and 

h +i = Q^f/c- 

So far, the process defined above is discrete in time. To obtain a time 
continuous process, we let the master vector V be a function of time, and the 
times between the jumps (collisions) exponentially distributed. In this way, 
if Fn(V, 0) is the probability distribution of the N particles on f) ;V 1 (x/N ) at 
time 0, the time evolution of Fn (V, t ) is described by the following master 
equation 

^dt'^ =^ f N(V,f), (3.3) 

where 

IC tp = N[Q^ - I] (3.4) 

and I is the identity operator. A more complete discussion of master equa¬ 
tions of this kind may be found in 0. 
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Note that, for <p(v) = v 1 , the collision operator K, ( p is the same as the 
collision operator K in the Kac model. We have propagation of chaos for 
the master equation (13.311 : 

Theorem 3.1. Assume that the family of initial densities {-Fn(W0)}ncn on 
n N ~\VN) is f(v,0)—chaotic. Then, the family of densities {Fn(T, OInciN/ 
where FN(V,t) is the solution to (I3.3P is f{v,t)—chaotic. Moreover, the density 
f(v, t ) satisfies the following equation 

r) r r 2n dO 

—f(v, t ) - [f(v(6), t)f{w(0), t ) - f(v, t)f(w, f)] (3.5) 

with f(v, 0) = fo(v ) and v(6), w(0) given by (I3.1D . 

Proof. The proof follows by the same arguments as in fgj. For a more de¬ 
tailed proof where propagation of chaos is shown for more general master 
equations, we refer to J2j|. □ 


4. Chaotic measures in higher dimensions 


In Section 2 we proved that the uniform density with respect to the mi- 
crocanonical measure on O N 1 ( y/N ) is Ce~ 2 °'W ! i -chaotic, v G IR. The goal 
of this section is to generalize this to the case v G IR 2 where now both the 
energy and momentum are conserved; the generalization to IR rf , d > 2 may 
be treated in the same way. The calculations here are formal. For p G IR 2 , 
define 


. N 

T n (Ve,p) = \ (v lr ...,v N ) G R 2N =2E, X> = p 

i =1 


N 

E 

i=l 


(4.1) 


We assume that E and p are chosen such that F' v (\/E, p) is non-empty. The 
classical case when <p(v) = |u| 2 has been thoroughly investigated in H3l . 

For p = (pi, pi ) and v, = (vn,Vii), i = 1,... N, a measure /iE,pi,p 2 con_ 
centrated on r N (\/E, p) is defined by 

N N N 

Pe, P lrP2 = 6{2E - J2f(0i))d(pi -^ l v a .)S(p 2 ~ ( 4 - 2 ) 

;=1 i =1 i =1 


The product of the Dirac measures is well defined since the hyper surfaces 
defined by setting the arguments of the Dirac measures to zero are mutually 
transversal. Let 


Z(E,p 1 ,p 2 ) =f 
Jr 2N 


N N N 

S(2E -J^(p(vi))S(p 1 - Y^ v n)d(P2 -Y^v a )dv i 

i=i i=i i=i 



As in Section 2, we need to determine the asymptotic behaviour of Z(E, p\, pi) 
with E = N for large N. We note that in the case of relativistic collisions, 
with cp(v ) = ^/\v \ 2 + l — 1, the measure p-E,m,v 2 is Lorentz invariant (see 
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In the sense of distributions, the 5 function is the inverse Fourier trans¬ 
form of the function 1 and formally can be written as 

5{x) = ^- [ e^d£. (4.4) 

271 J R 

For z = (zi, Z2, z 3 ), we can formally write 


Z(E, pi, p 2 ) 

= - -L, [ [ e ^2E-lZ 1 <f(^)) Z 3 e Kpi-E£i^)zi e 4P2-E l ! li^2) Z 2 rfz dz dz dv 

(2/r ) 3 Jr 2 N J r3 

where dV = tto| ... dv\>. With E = N, the last equality is 
Z(N,p 1/ p 2 ) = 

f e i Pl z 1+ ip 2 Z 2 ( e 2iz 3 r g-wvifo-ivuzi-ivuzifo dv \ d Zl dz?dz 3 . 
(2n) 3 Jr3 V J rz J 

Note that here p\ and p 2 are assumed to be independent of N. A natural 
and straightforward variation is to replace p\ by Np\ and p 2 by Np 2 . 


Let 


(4.5) 


< 7 ( Zl ,z 2/ z 3 ) = e ir>lZl+ip2Z2 , 

S(z I,z 2 ,z 3 ) = 2 iz 3 + log/ R2 e~ ^ (®1 ) z 3 -ronzi -W12Z2 dvndvn, 
where q : C 3 —> R, S : C 3 —> R. With these notations, we can write 

Z(N,pi,p 2 ) = , 1 , 3 [ q(z\, z 2 , z 3 )e NS ^ Zl,Z2 ' z ^dz\dz 2 dz 3 . (4.6) 

(27t)- 5 2r3 

The right hand side of the last equality is a saddle point integral in dimen¬ 
sion 3. The asymptotic behaviour of Z(N, p\, p 2 ) for large N is determined 
by the saddle points of S(zi,z 2 ,Z 3 ), i.e., points (zi,z 2 ,z 3 ) such that 

VS(z i,z 2 ,z 3 )=0. (4.7) 

Using (I4.5D . we need to solve the following system of equations: 

i f R2 (pfale-tyfate-^-^dvudvu 


2i — 


J r2 - iv nz 2 dv n dv 12 

i f R 2 v n e-^M^-imzi-m^dvudvu 
f R 2 e-^M^-wuzi-w 12 z 2 dv n dv 12 

i f R 2 vne-^^-^-^^dvudvu 


= 0 , 

= 0 , 

= 0 . 


(4.8) 

(4.9) 
(4.10) 


f^e-W^-i^i-^dvudw 

Since ([> is even, it follows that Zi,z 2 = 0 are the unique solutions to equa¬ 
tions (14.9b and (14.101) . We can now obtain z 3 by solving the following equa¬ 
tion: 

[ (2 — cb(vi))e~ tc P( l ’ 1 ' >Z3 dvi = 0. (4.11) 

Jr 2 
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Assuming that (0, 0 , 23 ) is the unique solution to (14.71) . and that we can de¬ 
form the integration domain in (14.61) to contain (0,0, Z3), we find that 


Z(N, pi, p 2 ) ~ 


1 / 27r\ 3/2 1 

(2 tt)3 VaTJ (detS"((0,0,z 3 ))) 1/2 

1 _I_ e 2Niz 3 | 

( 2 /rN ) 3 / 2 (detS"((0,0,z 3 ))) 1/2 


e NS( 0,0,23)^0,23) 

(4.12) 


By the discussion in Section 2, this procedure shows formally that the uni¬ 
form distributions on T N (VN,p) with respect to the measure PN,p lr p 2 is 
Ce~ z -chaotic, where Zq = iz 3 and £3 the unique solution to (14.111) and 


f R2 e z o<P( v )dv 

Sznitman's method, referred to in Section [T| and the end of Section [ 2 } is 
not restricted to the one-dimensional setting, but could formally be used 
also here: with N spatial dimensions, we would have d = N + 1, and 
h in Equations dl.131) and (11.141) would be a generalized four-momentum, 
h = (z?i, z>2, ...,VN,(p(v 1, ...,Vn)). However, the same difficulty as in the one¬ 
dimensional case would appear here, and a modification of Sznitman's ar¬ 
gument would be needed, in the same way. 


5 . Appendix 

The saddle point method is used to determine the asymptotic behaviour 
of integrals depending on a parameter. For a detailed description we refer 
to | 8 j. Without proof we only state below the saddle point method which is 
concerned with this paper. 

One-dimensional saddle point method 

Let 7 be a contour in the complex plane. Assume that q and S are ana¬ 
lytic functions in a neighborhood of the contour 7 . Consider the following 
integral 

F(A) - [ q(z)e xs ^dz. (5.1) 

■n 

A point zo G C is called a simple saddle point of the function S : C —> C if 
S'(z 0 ) = 0 and S"(zo) 7 ^ 0. Assume that zo G 7 is the unique simple saddle 
point of S. Then as A —>• 00 

F < A )=\/51i eASl2o) (^) +o <x0' (5 - 2) 

If there are more than one saddle point, F will be expressed as a sum over 
these points. 

Many-dimensional saddle point method 
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Let 7 be an N-dimensional smooth compact manifold. Consider the follow¬ 
ing integral 

F( A) — f (j(z)e AS ^ z - ) dz. (5.3) 

Jj 

where z = (z \,..., Z\i) G C ;V and the functions ij(z) and S(z) are assumed 
to be analytic in a domain D containing 7 . A point Zo is called a simple 
saddle point of S(z) if VS(zo) =0 and detS"(z 0 ) 7 ^ 0. Assume that Zo G 7 
is the unique simple saddle point of S. Then as A —t 00 

F < A > = (f) N/2 (de tS 4o)T/^ M (’<*> + ° ( I>) ' <5 ' 4) 

If there are more than one saddle point, F is be expressed as a sum over 
these points. 
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